DOC  FILE  COPY,  A[)  a 0 3 8 0 7 8 


AfOSW  - TR-  7 y 


7 

/ 

1 0 2 5 G 

FINAL  TECHNICAL  REPORT 
to  the 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

from 

Eugene  Herrin 
and 

Tom  Goforth 

Dallas  Geophysical  Laboratory 
Southern  Methodist  University 
Dallas,  Texas  75275 


Approved  for  r-  ' ' ' ' r ' " a : 

distribution  unliaittd. 


ARPA  Order:  2382 

Program  Code:  4F10 

Name  of  Contractor:  Southern  Methodist  University 

Effective  Date  of  Contract:  February  1,  1973 

Contract  Expiration  Date:  July  15,  1976 

Total  Amount  of  Contract  Dollars:  $711,731 

Contract  Number:  F44620-73-C-0044,  P00005 

Principal  Investigator  and  Phone  Number:  Eugene  Herrin 

AC  214  692-2760 

Program  Manager  and  Phone  Number:  Truman  Cook,  Director  of 

Research  Administration 
AC  214  692-2031 


Title  of  Work:  The  Detection  and  Research  of  Rayleigh  and 

Love  Waves 

University  Account  Number:  80-58 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  2382 


D D C 

Effsn^QiZE 


APR  8 1977 

JUElSEirD'E 


D 


1 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  ( AFSS) 
r •:  CE  OF  TRAN SK I T IAL\JO  DDC 

• technical  report  has  been  reviewed  and  is 
i.  - ovt’d  for  public  release  IA'.V  AFR  ISO-12  (7b)* 
D-lS trifcution  is  unlimited. 

A.  D.  BLOSE 

Technical  Information  Officer 


% 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fW7-»n  Dtf  Enl.red) 

! REPO?T  DOCUMENTATION  PAGE 


nr  ndm 


2 GOVT  ACCESSION  NO. 


AI-OSK:  - TR-  7 7"  V 2 5 {',  / 


«.  I Tl  Tl  F /Anri  SuhllllcJ 

I XHE  DETECTION  AND  RESEARCH  OF  RAYLEIGH  AND 
i LOVE "WAVES#  JX 


7. — Au  TMOR£a>  . 

! Eugene  Herrin  -V." 

Tom  Goforth 

— -•  / 

9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Dallas  Geophysical  Laboratory  / 

Southern  Methodist  University  ^ 

Dallas  TX  75275 

tl.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

ARPA 

1400  Wilson  Blvd. 

Arlington  VA22209 

14  MONITORING  AGENCY  NAME  A AODRESSrff  differin'  Iron,  Controlling  Ofllct) 


READ  instructions 
BEFORE  COMPLETING  FORM 

3 RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  t PERIOD  COVERED 


/ Final  , / . „ 


€ PERFORMING  ORGT’ 


AFOSR/NP 

Bolling  AFB , Bldg. #410 
Wash  DC  20332 


e CONTRACT  OR  GRANT  NUMBER^.) 

f F4462tf-73-C-0Q44,  

S-  : s: 

It)  PROGRAM  ETtrwENT  PROJECT,  TASK 
-AREA  A WORK  UNIT  NUMBERS 

&Z7uiE  4F1U 
2382-06 

tz  REPORT  O ATE  -1 

ij  I Feb  4977  ! * 

~13.  NUMBER  of  ^ ' 


15.  SECURITY  CLASS,  (ot  this  report) 


Unclassified 

15a.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


M6-  DISTRIBUTION  STATEMENT  (ot  thin  Report) 


Approved  for  public  release? 
distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (ot  the  ebstrect  entered  In  Block  20,  It  dltterent  from  Report) 


is  supplement  ary  notes 


TECH  OTHER 


I 19  KEY  WORDS  (Continue  on  reverse  side  It  neceesery  end  Identify  by  block  number) 


7J\ 


20  ABSTRACT  (Continue  on  reverse  side  It  necessery  end  Identify  by  block  number) 

‘'.Phase-matched  flitters  are  defined  as  a class  of  linear  filters  in  which  the 
Fourier  phase  of  the  filter  is  made  equal  to  that  of  a given  signal.  An 
Interative  technique  is  described  which  can  be  used  to  find  a phase-matched 
filter  for  a particular  seismic  signal.  The  process  is  then  applied  to  digital 
records  of  Rayleigh  waves  from  a synthetic  source  with  propagation  across  55 
deg.  of  continental  path,  an  earthquake  in  the  Greenland  Sea  recorded  in  Texas, 
and  a nuclear  explosion  in  Novaya  Zemlya  recorded  in  New  Mexico.  ^ h " >■  ' 


DD  ,^73  1473  EDITION  OF  I NOV  AS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dmtm  En(»r»rf) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  P AOEfUTin  Oda  gnfrtd) 


Application  of  the  filter  allows  multiple  arrivals  to  be  identified  and 
removed  and  allows  recovery  of  the  complex  spectrum  of  the  primary  wave 
train  along  with  its  apparent  group  velocity  dispersion  curve.  The 
amplitude  spectrum  of  the  primary  signal  obtained  by  this  linear  spectrum 
of  the  primary  signal  obtained  by  this  linear  process  is  not  contaminated 
by  interference  from  multipath  arrivals.  The  filtering  process  also 
provides  significant  improvement  in  signal-to-noise  ratio;  greater  than 
factor  of  four  for  the  Greenland  Sean  and  Novaya  Zemlya  events. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOC|IFA«n  Enf •red) 


BESmmmvm 


'f 


PHASE-MATCHED  FILTERS: 
Application  to  the  Study 
of  Rayleigh  Waves 


D D C 


I | (1L1 

U 


APR  8 1377 


uyibisEurnE 

D 


Eugene  Herrin  and  Tom  Goforth 
Geophysical  Laboratory 
Southern  Methodist  University 
February  1977 


I 

* 

* 


DISTRIBUTION  STATIC :LWT_A 

Approved  to i public  roieQ3c: 
Distribution  Unlkaitod 


S 

i 


ABSTRACT 


Phase-matched  filters  are  defined  as  a class  of  linear 
filters  in  which  the  Fourier  phase  of  the  filter  is  made 
equal  to  that  of  a given  signal.  An  iterative  technique  is 
described  which  can  be  used  to  find  a phase-matched  filter 
for  a particular  seismic  signal.  The  process  is  then  applied 
to  digital  records  of  Rayleigh  waves  from  a synthetic  source 
with  propagation  across  55  deg.  of  continental  path,  an  earth- 
quake in  the  Greenland  Sea  recorded  in  Texas,  and  a nuclear 
explosion  in  Novaya  Zemlya  recorded  in  New  Mexico.  Application 
of  the  filter  allows  multiple  arrivals  to  be  identified  and 
removed  and  allows  recovery  of  the  complex  spectrum  of  the 
primary  wave  train  along  with  its  apparent  group  velocity 
dispersion  curve.  The  amplitude  spectrum  of  the  primary 
signal  obtained  by  this  linear  process  is  not  contaminated 
by  interference  from  multipath  arrivals.  The  filtering 
process  also  provides  significant  improvement  in  signal-to- 
noise  ratio;  greater  than  a factor  of  four  for  the  Greenland 
Sea  and  Novaya  Zemlya  events . 
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INTRODUCTION 

Matched  Filters 

The  matched  filter  has  been  defined  by  Turin  (I960, 
p.  311)  as  follows:  Given  any  physical  waveform,  s(t),  then 

the  filter  matched  to  s(t)  is,  by  definition,  one  with  the 
impulse  response 

M-T)  = ks(A-r),  (1) 

where  k and  A.  are  arbitrary  constants.  In  equation  (1),  when 
we  consider  •£■  to  have  the  units  of  time,  then  A denotes  a 
shift  in  time.  We  can  consider  the  factor  k to  be  1 and  apply 
a shift  in  the  time  axis  so  that  A =0.  Thus  the  impulse  re- 
sponse of  the  filter  matched  to  s(t)  is  simply  s (-  f ).  Con- 
volution of  the  matched  filter  with  the  signal  is 

h(T  ) * s (t)  = s ( — T ) * s (t)  =p>js  (o)|  2 exp  i [-  <r  (4 >)+  «-(o)J  (2) 
where  the  symbol  means  "...  has  the  Fourier  transform  . . ." 
Thus  the  right  hand  side  is  just  jS(to)J  2.  Now  consider  the 
cross-correlation  of  s(t)  with  its  matched  filter  having  the 
direction  of  time  reversed 

s (t)  <$£)  s (t)  ■=^|s(<0)|2  exp  i j\r(  cj  ) - o~(cj  )j 
which  is  recognized  as  the  auto-correlation  of  s(t).  We  see 
that  the  output  obtained  by  filtering  s(t)  with  s(-'f)  is  the 
same  as  that  obtained  by  auto-correlating  the  signal,  s(t). 
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Turin  (I960,  p.  312)  shows  that,  in  the  presence  of  "white 
noise",  the  matched  filtering  process  maximizes  the  signal 
to  noise  power  ratio.  In  the  following  discussion  we  will 
assume  that  the  instrument  used  to  record  seismic  data,  the 
seismograph,  has  an  amplitude  response  that  is  approximately 
the  inverse  of  the  background  noise,  and  thus  the  requirement 
for  "white  noise"  is  essentially  met. 

Phase-Matched  Filters 

We  consider  convolution  and  cross-correlation  of  a signal, 
s(t),  with  a time  function,  f(t),  to  be  represented  by  the 
following  operations  and  their  Fourier  transforms, 

s(  t)  * f(D^|s(«)l  |F(o)|  exp  i Jjs-  ( )j 

and 

s (t)  & f (t)  =t>|s(co)]  |f(63  ]|  exp  i[<s-(«)  (3) 

Now  suppose  that  we  choose  fp(t)  such  that  the  Fourier  phase  is 
the  same  as  that  of  s(t).  We  define  the  class  of  linear  oper- 
ators, fp(t),  such  that  <s-(co)  = <£^(60),  as  phase-matched 
filters  with  respect  to  the  signal,  s(t).  The  output  of  the 
second  of  the  above  operations  will  then  have  the  Fourier 
transform,  |s(&J)|  jFp(<5J  )|  , and  will  be  an  even  function  in 
the  time  domain  as  is  the  auto-correlation  function.  The  out- 
put of  the  first  operation,  convolution,  with  the  direction  of 
time  reversed,  that  is,  ) substituted  for  fp(f  ),  will 


3 


have  the  same  Fourier  transform.  We  call  this  output,  whether 
obtained  by  filtering  or  correlation,  a pseudo-autocorrelation 
function  (PAF) . If  means  can  be  found  for  matching  the  Fourier 
phase  of  the  signal  and  the  filter,  then  the  time  domain  out- 
put, the  PAF,  will  depend,  for  a given  signal,  only  upon  the 
amplitude  spectrum  of  the  particular  phase-matched  filter 
used  in  the  operation. 

Now  consider  some  possible  choices  forjFp(<w)|. 

|Fp(o)|  = |S(4))|  (a) 

|F  p(4>)|  = 1 (b) 

1 

]Fp(<o)|  = |S(A>))  (C) 

If  choice  (a)  is  made,  the  phase-matched  filter  becomes  the 
matched  filter  and  maximizes  the  signal  to  noise  power  ratio 
assuming  "white  noise".  If  we  were  using  a power  detector  to 
indicate  the  presence  of  a low-level  signal,  choice  (a)  should 
give  the  best  results.  In  this  case  the  PAF  becomes  the  auto- 
correlation function.  On  the  other  hand,  choice  (c)  is  simply 
the  inverse  filter  and  the  PAF  becomes  the  impulse  function. 

In  practice,  choice  (c)  would  maximize  the  time  resolution 
of  the  output,  but  would  greatly  reduce  the  signal  to  noise 
ratio  (Turin,  1960,  p.  318) . If  our  problem  is  one  of  low 
signal  to  noise  ratio  we  should  use  a filter  spectrum  approxi- 
mating choice  (a),  but  if  we  have  interfering  signals,  all 
with  large  signal  to  noise  ratios,  a spectrum  approaching 
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that  in  choice  (c)  would  be  more  appropriate.  In  the  following 
examples  we  have  used  choice  (b),  a convenient  compromise 
between  maximum  signal  to  noise  improvement  and  maximum  time 
resolution. 

The  above  relationships  result  from  the  use  of  infinite 
integrals  and  continuous  functions  in  the  time  and  frequency- 
domains.  The  results,  however,  hold  for  practical  applications 
involving  finite  sums  and  sampled  functions  provided  proper 
precautions  are  taken  to  avoid  wrap-around  and  aliasing  in 
computing  finite  Fourier  transforms. 


» 
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APPLICATIONS  TO  RAYLEIGH  WAVES 


Alexander  and  Lambert  (1971)  discuss  the  application  of 


matched  filters  to  the  analysis  of  surface  waves  and  to  the 


determination  of  their  amplitude  spectra.  They  show  that  by 


"pre-whitening"  the  matched  filter  the  length  of  the  output 


function  was  reduced  by  about  a factor  of  two.  They  state 


(p.  7)  that  "the  only  negative  feature  of  the  whitened 


matched  filter  is  that  as  a detector  it  will  be  less  satisfac- 


tory than  the  unwhitened  filter."  In  fact,  they  were  using 


a phase-matched  filter,  as  defined  in  the  previous  section. 


with  an  amplitude  spectrum,  (F(CJ)j  , approximately  equal  to 


1 over  the  band  of  interest.  The  length  of  the  PAF  was  halved 


thus  the  time  resolution  was  twice  that  of  the  matched-filter. 


Alexander  and  Lambert  (1971)  did  not  show  how  to  obtain 


the  matched  filter  when  the  exact  form  of  the  signal  is  un- 


known. In  this  section  we  describe  an  iterative  technique 


which  can  be  used  to  find  a phase-matched  filter  for  a given 


seismic  signal.  The  process  is  then  applied  to  digital  records 


of  Rayleigh  waves. 


Table  1 is  a flow  diagram  of  the  computer  code  we  use  in 


the  phase-matching  process.  We  assume  that  the  seismic  signals 


result  from  a dispersive  process  with  a continuous  dispersion 


TABLE  i 
FLOW  DIAGRAM 


FREQUENCY  DOMAIN 
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curve  (group  velocity  vs.  frequency).  The  group  delay,  tgr ( a ) , 
associated  with  the  signal  is  the  epicentral  distance  divided 
by  the  group  velocity  at  frequency  6)  . From  Papoulis  (1962, 
p.  134)  we  note  that  the  group  delay  and  Fourier  phase  of  the 
signal  are  related  as  follows: 

d 0(d>) 


tgr(6>) 

= 

d 

(4a) 

0 (4ix) 

= 

1 tgr  ( &) ) d 0$  . 

and 

(4b) 

In  the  code  the  signal  is  Fourier  transformed  to  obtain 
the  amplitude  and  phase  spectra.  These  finite  transforms  in- 
volve only  positive  times  and  frequencies  so  that,  in  practice, 

T r T 1 

time  shifts  of  ^ and  phase  shifts  of  exp  1 -i  ^ O , w^ere  T is 
the  length  of  the  time  domain  window,  must  be  taken  into  account. 
We  provide  as  input  a trial  group  velocity  dispersion  curve 
and  an  amplitude  spectrum  for  the  filter.  Using  the  epicentral 
distance  to  the  signal  source  and  Eq.  (4b)  we  compute  the 
Fourier  phase  of  f(t)  and  perform  the  correlation  in  the  fre- 
quency domain  (Eq.  3).  The  result  is  then  transformed  to  the 
time  domain,  windowed  to  reject  correlation  functions  from  inter- 
fering signals  or  multipath  arrivals,  then  transformed  again. 

The  result  will  have  the  complex  spectrum 

|s(d>)|  |f  ( &)  ) J exp  i [<T(<0)  - 4>(00)J  ' (5) 

We  then  "unwind"  the  difference  in  Fourier  phases,  <s~  (W  ) - <p  (CJ  ) , 
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which  in  practice  is  multiplied  by  the  linear  phase  shift 
associated  with  the  offset  in  the  lag  domain  window,  and 
obtain  a smoothed  estimate  of  its  derivative.  The  smoothing 
process  assumes  that  the  first  derivative  of  the  dispersion 
curve  is  continuous;  hence,  by  Eq.  4a,  the  second  derivative 
of  the  Fourier  phase  spectrum  must  be  continuous.  This  esti- 
mate is  used  to  correct  the  group  delay  of  the  trial  filter 
(Eq.  4a) . The  process  is  repeated  until  the  phase  spectra  of 
the  filter  and  the  desired  signal  in  the  band  of  interest  are 
identical.  We  have  then  obtained  the  phase-matched  filter, 
fp(t),  having  the  given  amplitude  spectrum,  jFp(Cj)J  . The 
spectrum  of  the  signal  can  be  recovered  because  the  Fourier 
transform  of  the  resulting  PAF  is  J S ( CO ) J Jf^ ( 6J ) . In  the 
following  examples  the  filter  spectrum  is  "whitened";  that  is, 
jFp(co  )j=  1 in  the  frequency  band  of  interest  (periods  of  10  to 
100  seconds). 


and  3 represent  reflections.  The  source-strength  is  assumed 
to  vary  azimuthally  so  that  arrivals  1 and  2 have  comparable 
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strength.  Arrival  3 has  about  *5  the  amplitude  of  the  other 
two.  The  time  separations  are  100  sec  between  arrivals  1 
and  2 and  300  sec  between  1 and  3.  The  same  dispersion 
curve  was  used  to  generate  all  three  arrivals,  but  because 
of  the  difference  in  path-lengths  the  group  delays  and 
Fourier  phases  of  the  arrivals  will  be  different.  Fig.  2 
shows  the  synthetic  signal  for  an  epicentral  distance  of  55 
deg.  The  Airy  phase  for  the  second  arrival  can  be  seen  but 
the  end  of  the  third  arrival  is  not  shown. 

The  trial  dispersion  curve  used  as  input  to  the  phase- 
matching code  was  significantly  different  from  the  one  used 
to  generate  the  synthetic  signal.  Fig.  3 shows  the  correlation 
functions  generated  in  four  iterations  of  the  code.  An  iter- 
ative procedure  is  required  because  the  smoothing  process 
used  in  the  estimation  of  the  Fourier  phase  derivative  is 
numerically  damped  and  because  the  operator  may  change  control 
parameters  during  the  calculation.  Pass  1 in  Fig.  3 shows 
the  correlation  function  obtained  using  the  trial  dispersion 
curve.  Multipathing  is  definitely  indicated  so  the  operator 
chose  a fairly  narrow  lag  window  centered  on  the  predicted 
arrival  time  of  Rayleigh  energy  with  greater  than  100  sec. 
period.  The  window  was  narrowed  on  subsequent  passes,  and 
in  Pass  4 a nearly  perfect  PAF  was  obtained;  that  is,  the 
phase  of  the  primary  signal  was  matched  within  the  limitations 
of  accuracy  imposed  by  the  sampling  rate  (1  per  sec.).  After 
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the  final  iteration  essentially  all  the  primary  energy  has 
been  compressed  into  about  50  sec.  of  lag  time.  The  correla- 
tion functions  for  arrivals  2 and  3 show  that  the  filter  is 
not  phase-matched  to  these  arrivals  because  their  group  delays 
are  different  from  that  of  the  primary.  In  Fig.  3,  Pass  4, 
note  the  correlation  functions  for  the  second  and  third  arri- 
vals are  displaced  to  the  right  of  the  actual  arrival  times, 
as  shown  by  the  dashed  lines,  because  of  the  longer  travel 
paths.  This  effect  helps  in  extracting  the  primary  energy 
when  the  separation  times  between  arrivals  becomes  less  than 
about  50  sec.  The  amplitude  spectrum  of  the  primary  was 
recovered  by  taking  the  Fourier  transform  of  the  windowed 
correlation  function  of  Pass  4. 

Fig.  4 shows  the  starting  (solid  line)  and  final  (dashed 
line)  dispersion  curves  for  the  primary  arrival  in  the  synthetic 
signal.  The  dashed  curve  agrees  with  the  one  used  to  construct 
the  synthetic  signal  within  +0.005  km/sec.  This  test  demon- 
strates the  use  of  the  filtering  process  to  resolve  overlapping 
signals,  but  without  the  presence  of  noise.  The  next  example 
uses  data  from  an  earthquake  with  a significant  amount  of  back- 
ground noise  present. 

Greenland  Sea  Event 

This  particular  signal,  recorded  from  a vertical  seismometer. 


GROUP 

VELOCITY 

km/sec 


3.0 


2.0 


1.0 


i 1 ^ 

10. 

PERIOD 


-t 

100. 


-I — I — t- 


Fig.  4 


100C 


The  starting  (solid  line)  and  final  (dashed  line)  dis 
persion  curves  for  the  synthetic  signal. 
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was  chosen  because  it  shows  very  regular  dispersion  with 
almost  no  indication  of  multipath  interference.  Table  2 
gives  the  pertinent  data  on  this  earthquake.  The  record 
was  made  in  a salt  mine  at  Grand  Saline,  Texas,  using  a 
long  period  seismometer  with  an  amplitude  response  that  is 
essentially  the  inverse  of  the  background  noise  from  periods 
of  10  to  100  sec.  The  spectra  shown  in  this  paper  have  not 
been  corrected  for  instrument  response. 

A trial  dispersion  curve  was  selected  and  an  initial 
filter  was  calculated  in  the  frequency  domain.  This  filter 
has  an  amplitude  spectrum  that  is  flat  from  periods  of  4 sec. 
to  100  sec.,  thus  the  variations  in  amplitude  seen  in  the  time 
domain  representation  of  the  filter  result  from  the  properties 
of  the  dispersion  curve  (see  Fig.  6) . 

Fig.  7 shows  the  correlation  functions  obtained  in  four 
iterations  of  the  code.  As  the  group  delay,  and  hence  the 
dispersion  curve,  was  modified  the  correlation  function 
approached  a perfect  PAF  and  an  excellent  phase-match  was 
obtained.  Fig.  8 shows  the  starting  (solid  line)  and  final 
(dashed  line)  dispersion  curves  for  the  Greenland  Sea  event. 

It  must  be  emphasized  that  the  final  curve  is  an  apparent, 
group  velocity  dispersion  curve.  We  have  not  corrected  the 
curve  for  the  phase  response  of  the  instrument;  furthermore, 
our  studies  on  data  from  long  period  arrays  (LASA  and  ALPA) 
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Table  2 

Data  pertinent  to  the  earthquake 
in  the  Greenland  Sea 


ORIGIN  TIME 
LATITUDE 
LONGITUDE 
DEPTH 

MAGNITUDE  (ir^) 
RECORDING  SITE 
EPICENTRAL  DISTANCE 


1970-139-02-07-41 . 5 
N 79.2  Deg. 

E 2.5  Deg. 

N 

4.8 

Grand  Saline,  Texas 
59.70  Deg. 


Table  3 

Data  pertinent  to  the  nuclear  explosion 
in  Novaya  Zemlya,  U.S.S.R. 


ORIGIN  TIME 
LATITUDE 
LONGITUDE 
DEPTH 

MAGNITUDE  (1%) 
RECORDING  SITE 
EPICENTRAL  DISTANCE 


9 

I 

* 


s 


1975-294-11-59-57.3 
N 73.35  Deg. 

E 55.08  Deg. 

0 

6.5 

Albuquerque,  New  Mexico 
71.28  Deg. 


The  vertical  component  of  the  Rayleigh  wave  recorded 
from  the  Greenland  Sea  event  (see  Table  2) . 
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The  starting  (solid  line)  and  final  (dashed  line)  dis 
persion  curves  for  the  Greenland  Sea  event. 
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show  that  even  the  primary  arrival  of  Rayleigh  energy  may 
not  follow  a great  circle  path.  These  problems  will  be 
treated  in  detail  in  subseguent  papers. 

In  Fig.  9 the  amplitude  spectrum  of  the  windowed 
correlation  function  from  Pass  4 (solid  line)  is  superimposed 
on  the  spectrum  (dashed  line)  obtained  by  transforming  the 
entire  signal  shown  in  Fig.  5.  There  is  some  smoothing  of 
the  latter  spectrum  caused  by  the  narrow  lag  window  but  the 
spectra  are  essentially  the  same.  The  spectral  hole  at  a 
period  of  about  20  sec.  is  real;  that  is,  the  hole  was  not 
caused  by  multipath  interference  and  is  presumably  indicative 
of  the  source  spectrum.  The  filtering  process  has  resulted 
in  signal-to-noise  improvement  of  greater  than  a factor  of 
4 in  the  band  of  interest  (compare  Figs.  5 and  7). 

Novava  Zemlva  Event 

Fig.  10  shows  the  signal  from  the  vertical  module  of 
a downhole,  long  period  seismograph  at  Albuquerque,  New  Mexico 
(U.S.G.S.),  for  Rayleigh  waves  from  an  explosion  in  Novaya 
Zemlya,  U.S.S.R.  Pertinent  data  for  this  event  are  given  in 
Table  3.  The  signal  shows  a beat  pattern  and  an  apparent 
repetition  of  the  Airy  phase  which  are  indicative  of  multipath 
interference.  A trial  dispersion  curve  for  the  north  polar 
path  was  chosen,  and  a whitened  amplitude  spectrum  was  then 


'l.o  JO.  □ 100.0  1000.0 


SEC 


The  signal  spectrum  (dashed  line)  and  the  spectrum  of 
the  final  correlation  function  (solid  line)  for  the 
Greenland  Sea  event. 


The  vertical  component  of  the  Rayleigh  wave  recorded 
from  the  Novaya  Zemlya  event  (see  Table  3). 
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used  to  construct  the  filter  and  perforin  the  first  correlation. 
The  resulting  correlation  function  and  those  for  two  subsequent 
passes  are  shown  in  Fig.  11.  By  Pass  3 an  excellent  PAF  was 
obtained.  The  starting  (solid  line)  and  final  (dashed  line) 
dispersion  curves  are  shown  in  Fig.  12. 

There  is  a good  indication  of  a second  arrival  80  to  90 
sec  after  the  primary  and  another  arrival  about  170  sec  later. 

Of  course,  the  filter  was  not  phase-matched  to  these  later 
arrivals.  The  PAF  from  Pass  3 was  windowed  in  the  lag  domain 
so  as  to  exclude  the  second  arrival  but  retain  all  of  the  primary 
correlation  function.  The  signal  spectrum  (dashed  line)  and 
amplitude  spectrum  of  the  windowed  PAF  (solid  line)  are  shown 
in  Fig.  13.  The  ripple  in  the  signal  spectrum  caused  by  multi- 
path  interference  has  been  removed  by  the  filtering  process. 


Note  that  in  Fig.  11,  Pass  3,  the  total  width  of  the 
PAF  is  only  about  60  sec.  so  that  a second  arrival  only  30 
to  40  sec.  after  the  primary  could  have  been  removed  using 
the  final  filter.  Even  greater  resolution  may  be  obtained 
by  modifying  the  amplitude  spectrum  of  the  filter  so  that 

. 

the  PAF  approaches  the  band-limited  impulse  function  as 
,■*  discussed  in  the  first  section  of  this  paper. 
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DISCUSSION 


The  removal  of  the  effects  of  multipaths  using  phase- 
matched  filtering  provides  Rayleigh  wave  spectra  which  are 
much  more  closely  related  to  source  characteristics  than 
are  those  obtained  from  Fourier  transforms  of  the  entire 
signal,  and  at  the  same  time  results  in  a substantial  improve- 
ment in  the  signal  to  noise  ratio.  This  process  should  pro- 
vide the  data  necessary  for  studies  relating  surface  wave 
spectra  to  focal  parameters;  give  better  estimates  of  surface 
wave  magnitudes  and  their  relationship  to  explosion  yield; 
and  permit  re-evaluation  of  marginally  successful  techniques 
for  explosion/earthquake  discrimination  such  as  depth  deter- 
mination from  spectral  nulls  and  the  ratio  of  spectral  ampli- 
tudes at  periods  of  20  and  40  seconds. 

The  whitened  filter  used  on  the  Greenland  Sea  and  Novaya 
Zemlya  events  is  but  one  in  a continuum  of  possible  phase- 
matched  filters.  This  filter  appears  to  be  a useful  compromise 
between  matched  and  inverse ' filters  when  the  signal-to-noise 
ratio  is  moderate  and  the  secondary  arrivals  lag  the  primary 
by  at  least  30  to  40  seconds.  It  has  the  obvious  advantage 
of  requiring  no  prior  knowledge  of  the  true  amplitude  spectrum 


of  the  signal. 
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The  presence  of  multipaths  in  a surface  wave  train  will 
bias  conventional  methods  for  determining  apparent  group 
velocity  dispersion  curves  such  that  the  estimates  will  be 
too  slow.  The  removal  of  multipaths  by  phase-matched  filter- 
ing removes  this  bias  and  thus  improves  the  accuracy  of  the 
estimate.  The  precision  of  the  phase-matched  estimate  is  also 
improved  in  that  known  velocities  of  simulated  signals  can  be 
recovered  to  within  +0.005  km/sec. 

The  apparent  dispersion  curve  obtained  for  the  Greenland 
Sea  event  probably  represents  propagation  along  very  nearly 
a great  circle  path.  The  source-station  path  for  this  event 
consists  primarily  of  high-velocity  shield  and  mid-continent 
regions,  and  it  seems  likely  the  great  circle  and  least  time 
paths  coincide;  however,  there  is  no  reason  to  believe  that 
this  is  always  true.  Figure  14  shows  the  two  paths  and  indicates 
that  the  Novaya  Zemlya  - Albuquerque  great  circle  path  traverses 
a part  of  the  southern  Rocky  Mountains.  When  the  great  circle 
paths  lie  along  or  adjacent  to  major  structural  boundaries 
such  as  continental  margins  or  mountain  ranges,  the  possibility 
of  the  least-time  path  lying  along  another  route  must  be  con- 
sidered. Such  a circumstance  may  occur  for  many  source-station 
pairs,  and  obviously  complicates  the  problem  of  inverting 
apparent  dispersion  curves  to  obtain  earth  structure  since  the 
vertical  velocity  distribution  would  be  confounded  by  lateral 


A map  showing  me  = 
Sea  and  Novaya  Zemlya 
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inhomogenieties . Phase-matched  filtering,  with  its  ability 
to  resolve  overlapping  dispersed  wave  trains  into  a series 
of  impulse-like  functions,  offers  a new  approach  to  the 
surface  wave  inversion  problem.  A series  of  filtered  seis- 
mograms representing  various  paths  through  a region  can  be 
interpreted  in  terms  of  positions  of  lateral  boundaries  in 
much  the  same  way  that  deconvolved  seismic  reflection  and 
refraction  records  can  be  used  to  locate  velocity  interfaces 
at  depth.  The  complex  spectra  of  the  resolved  multiple 
arrivals  thus  obtained  can  then  be  related  to  changes  in 
physical  properties  at  the  interfaces. 

Phase-matched  filters  have  been  applied  to  Rayleigh 
waves  in  this  paper,  but  they  would  seem  to  have  general 
applicability.  Love  wave  signals  are  complicated  by  multi- 
path  Love  arrivals  as  well  as  off-axis  Rayleigh  waves,  hence 
phase-matched  filtering  should  prove  useful  in  their  study. 
The  application  of  these  filters  is  not  limited  to  signals 
which  are  strongly  dispersed  by  the  medium  through  which  they 
travel,  so  that  applications  to  the  study  of  body  waves  are 
anticipated. 
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